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The antiferromagnetic Heisenberg spin systems on the three-leg ladder are investi- 
gated. Periodic boundary condition is imposed in the rung direction. The system 
has an excitation gap for all antiferromagnetic inter-chain coupling (J_l > 0). The 
estimated gap for the strong coupling limit {J±/Ji oo) is 0.28Ji. Although the 
interaction is homogeneous and only nearest-neighbor, the ground states of the sys- 
tem are dimerized and break the translational symmetry in the thermodynamic limit. 
Introducing the next-nearest neighbor coupling (J2), we can see that the system is 
solved exactly. The ground state wave function is completely dimer-ordered. Using 
density matrix renomalization group algorithm, we show numerically that the original 
model ( J2 = 0) has the same nature with the exactly solvable model. The ground state 
properties of the ladder with a higher odd number of legs are also discussed. 
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§1. Introduction 

Nowadays there is much attention for quantum ladder systems.S^ One of the 
reason is that the ladder systems are the first step from one-dimensional toward 
two-dimensional.i'i'i-* The other reason is caused by the advance of experimental 
techniques. Many materials which seem to be ladders are found and investigated 
extensively.iiB The antiferromagnetic spin systems are among the simplest many 
body systems which include strong quantum effect. From the analysis using the 
0(3) non-linear a model,& antiferromagnetic spin ladders with an even number of 
legs are expected to have an excitation gap. In this CclSG, clS for the existence of 
an excitation gap, the system is not sensitive to the boundary condition perpen- 
dicular to the chain. For the two- leg S = 1/2 spin ladder there are many works 
and it is confirmed that there is an excitation gap.i'00) On the other hand, in the 
case of a ladder with an odd number of legs, the topology is important .0-' (For the 
fermion systems, there are works about the importance of the topology,Blli-' which 
treat the repulsive Hubbard ladder with weak electron-electron coupling limit us- 
ing field theoretical argument.) The odd-leg ladders with open boundary condition 
in the rung direction (we will call it as open ladder) has no frustrated interaction. 
These can be considered effectively as S* = 1/2 antiferromagnetic Heisenberg chain 
(AFHC)^©. On the other hand, periodic boundary condition in the rung direction 
(cylindrical ladder) causes frustrated interaction. The frustration changes the situ- 
ation dramatically. There is also attention for the dimer-ordered phase. We know 
the = 1/2 AFHC is unstable for the dimer interaction. The model becomes gapful 
as the dimer interaction is inserted.lii^ However the model breaks the translational 
invariance a priori. The Majumdar-Ghosh model is one of the typical models with 
which ground states are frustration induced dimer-ordered and have translational 
invariance in the finite size system.lli^ In this paper we show another example of the 
model with frustration induced dimer-ordered ground states. We consider following 
antiferromagnetic Heisenberg Hamiltonian, 

N' . . 

T^l = Jl^Yl ^n'n^ + ^J- X] 53 (1) 

71=1 i n=l i 

where H^^^ = Sn,i ■ Sm,j and Sn,i is S* = 1/2 SU{2) spin operator. We concentrate 
on the condition, Ji,J± > 0. Subscripts and superscripts n,m and i,j represent 
the site number in the rung direction and the chain direction respectively. The 
number of legs and the number of sites along the chain are and L respectively. 
(The total number of sites is NL.) For the open ladder we set A^' = — 1 and 
for the cylindlical ladder we set N' = N and A^ + 1 = 1. The boundary condition 
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along the chain direction is appropriately treated and denoted periodic boundary 
condition and open one as PBC and OBC respectively. In almost all part of this 
paper (except for §4), we will consider the three-leg ladder = 3. We add the 
next-nearest neighbor interaction, 

N' 

^2 = J2EE^^,t^^ (2) 

n=l i 

to see the system more explicitly. (See below.) Then the total Hamiltonian is 
7i = Til + 7^2- In the case of J± = 0, we know the model well for a few parameter 
points. The model with (Ji, J±, J2) = (-^1,0,0), which is exactly solved by Bethe 
An«atz,0 has no frustrated i,>teract,o,> a,>d gro„,>d state rs proved to be spm smglet 
by Marshall-Lieb-Mattis Theorem. The model has gapless excitation and belongs 
to the universality class of the level-1 SU{2) WZW model in the low energy limit. 
We know the long-distance asymptotic behavior of correlation functions, for exam- 
ple, spin-spin correlation function {Sq ■ Sr) and dimer-dimer correlation function 
{{So ■ Si){Sr ■ Sr+i)) — {Sq ■ Si){Sr ■ Sr+i) dccay as {—ly/r (up to logarithmic 
correction jli' ) . For the other parameter point (Ji,0, Ji/2), which is the so-called 
Majumder-Ghosh model,lli^ the ground states are well-known in the sense that we 
exactly know the wave function. The ground state wave function is dimerized and 
written down by direct products of two-site spin-singlet. Therefore the ground states 
are doubly degenerate and break the translational invariance. The existence of the 
excitation gap is exactly proved.© By the bosonization approach, we know that the 
perturbation by J2 is marginal. (Whether the operator is relevant or irrelevant is 
determined by the sign of the initial coupling constant. In the decoupled chains case 
(J_L = 0), the transition point is J2 = 0.2411 ■ ■ -.0^) On the other hand, perturba- 
tion by J_L is strongly relevant compared with that of J2-I^'l^ Therefore we expect 
that the total Hamiltonian with J_|_ > is essentially represented by the effective 
Hamiltonian with J±/Ji 00. 

Our strategy to investigate the model is as follows. First, we introduce the effective 
Hamiltonian in the strong coupling limit {Jj_ 00), which is the fixed point of the 
renomalizaion group. We add the next-nearest neighbor coupling ( J2) to the model. 
Then we show the model with special coupling (J2/J1 = 1/2, J_l/Ji — 00) is solved 
exactly. (In the sense that we can construct the exact ground state wave function.) 
Using the density matrix renomalization group (DMRG) algorithm introduced by 
White,!!-' we see the original model (J2 = and < J_l/Ji < 00) has the same 
nature with the exactly solved model by calculating the expectation value of the 
local Hamiltonian operator. We also confirm that there is no transition between 
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< J2/J1 ^ 1/2 through seeing no level crossing in low-energy excitation spectra 
by numerical exact diagonalization. The estimation of the excitation gap for various 
value of is also showed. In §4 we shortly mention the cylindrical ladders with a 
higher odd number of legs. For general odd number of legs A^, we construct effective 
Hamiltonian in the strong coupling limit, which has the same form with the case of 
AT = 3. 

§2. Strong-Coupling Effective Hamiltonian 

We first consider the strong coupling limit ( J_l 00) of the three-leg ladder. The 
reason is as following. In bosonization language (weak coupling analysis), the inter- 
chain coupling causes the relevant operator through the staggered part of the spin 
operator, 00) so generally the system is expected to flow into the strong coupling 
limit in renormalizaion processes, i.e. the strong coupling limit is the flxed point 
of the renomalization flow. In the case of the open ladder with Ji = (which is 
understood as Jj_ = 00), the ground state of each rung is a doublet. Then the 
effective Hamiltonian of the original Hamiltonian (|l|) in the strong coupling limit is 
represented by, 

'HeS = Jl'^^Si ■ S'j+1. (3) 

i 

Therefore, in this limit, the system can be described by 5 = 1/2 AFHC as noted in 
the introduction. The model has no excitation gap, power-law decaying spin-spin 
correlation, and it is classifled into the same universality of level-1 SU{2) WZW 
model.ii) The vanishing excitation gap is proved by Affleck and Lieb.B On the 
other hand, in the case of the cylindrical ladder, we cannot prove the existence of a 
gapless excitation by the theorem because the uniqueness of the ground state is not 
proved. The ground state of the three-site ring is four-fold degenerate, these states 
are. 
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where u = exp(^). The indices L,R represent the momentum of the three site 
ring, k = 27r/3, and — 27r/3 respectively. The four- fold degeneracy essentially dis- 
tinguishes the cylindrical ladder from the open ladder. The other states have higher 
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energy of order J^. Therefore in the strong couphng hmit J_l/ Ji oo, the effective 
Hamiltonian is 

Kfr = 4E^eff^'(«) with « = 4. (4) 

where Hlg'^^{a) = Si ■ Si+i (l + ^(rj^rj^i + Tj^Tj^i)). Though we write a exphcitly 
for convenience through this section, it should be considered as 4. The operators 
exchange indices L,R such as, 

r+|-L) = 0, T-\-L) = \-R), 

T+\-R) = \-L), T-\-R) = 0. 

(In these equations the dots represent | or t.) We define the diagonal operator 
which, with r^' = (r+ + r~)/2 and = (r+ — r^)/2i, satisfies commutation relation 
of three generators of SU{2). 
For the two-site system, the ground state wave function of the Hamiltonian (^) is, 

= (I Ti) - I it)) • (1^^) + 1^^)), (5) 

the corresponding energy is — 1(1 + a)Ji. In the following, we write this wave 
function as, •->-•. The wave function can be also written down in the three- leg 
ladder language. The wave function in the three-leg ladder is the sum of three terms 
which are [li, 2i] (8> [I2, 13] ^ [22, 23] and ones translated in the rung direction i.e.. 
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Here [■,■] represents the two-site spin-singlet state (| ti) ~ I it))- Of course the 
expression is not unique, we are able to write it in more complicated form with two 
site spin singlet states. 

Next we add the next-nearest neighbor coupling (Hamiltonian 7^2)- Then the total 
effective Hamiltonian is, 

= f E H^^'ic^ + f E H^'ic^- (6) 
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With special coupling Ji = 2J2, we find exact ground states as follows. Define the 
operator P, 

^..+i.+2 = l{Hi;^\a) + Hlt'^^^\a) + H%^\a)) + ^(1 + a). (7) 
Using this operator P the Hamiltonian is represented by 

n.s = ^Y.iP^,^+w-li^ + o^))- (8) 

The operator Pj j+i ,(_(_2 is small dimensional, we, therefore, easily proved the semi- 
positive definiteness of it. The number of eigenvectors with the eigenvalue zero is 
12, and these eigenvectors are represented in the following three types. 
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i.e. P is the linear combination of projection operators with positive coefficients, 
which projects out these three states. The situation is closely similar to the case of 
the Majumdar-Ghosh model. The difference is that P is not represented by a single 
projection operator. Then the two states, 

1^1)= . . ^ ^. ^ . . 

1 L 

are, at least, among the ground states with energy — 1(1 + a)LJi. ( For PBC case, 
the site 1 and L of the wave function \il)2) are connected by a line. For OBC case, \il)2) 
is degenerate because of the degrees of freedom of site 1 and L.) These ground states 
are evidently dimerized. In finite-size system with PBC we construct two ground 
states with the momentum and vr as + |'?/'2))/"\/2 and — \il)2))/\/^- We 
also confirmed that only these two states are ground states by numerical exact diag- 
onarization up to 12-site system. For the two ground states, obviously, correlation 
functions {SqS^) and (tqT^) have zero correlation except for the nearest neighbor. 
The wave function is site local, it strongly suggests the existence of the excitation 
gap. We, unfortunately, have not found the exact proof of it. 
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In the next section we perform the numerical calculation using DMRG algorithm. 
The DMRG prefers OBC to PBC Therefore we prepare the appropriate order pa- 
rameter for the discussion of the above models, which should show whether the 
ground state of the system is dimerized. Note that for the case of OBC, if we take 
the Hamitonian as, 

i=l " i=l 

= y E (nm.+2 - ^(1 + «)) + ^{H'Jia) + H^^''\a)l (9) 

the ground state is not degenerate for the finite-size system because of the last 
braketted terms. \1lJ2) has higher energy than \ipi) with ^(1 + a)Ji for all finite 
system size. Define the local Hamiltonian operator, 

= Hi:t\a) - -1- Y.{H^^\a)). (10) 

In the three-leg ladder case, we replace H^^^^^ with X]n=i(-f^n + -f^n,n+i + -^n^+*^^)- 
For the Hamiltonian (^), the expectation value of this operator apparently has long- 
range order, i.e., 

(^^;ti) = (-ir^(i+«), 

where we assume L = 2 mod 4. Therefore non-vanishing expectation value of 
this operator toward the center of the system explicitly means that the system has 
dimerized ground state, and breaks the translational symmetry in the thermody- 
namic limit. On the other hand, in the case of the open ladder (effectively this 
model is considered as S = 1/2 AFHC with OBC), the expectation value of this 
operator is decaying as (— 1)*"/ ^/r (up to logarithmic correction), where r is distance 
from the boundary. This is because the scaling dimension of the dominant part of 
the local Hamiltonian is one-half. 

§3. Numerical Calculations 

In this section we perform numerical calculations using DMRG algorithm to show 
the original three-leg ladder system has the same nature with the exactly solvable 
model, such as a non-zero excitation gap, a finite correlation length of the spin-spin 
correlation function and a doubled unit cell of the ground state wave function. In the 
following calculations, we use the finite-size algorithm of DMRG which is necessary 
for the system with a large number of degrees of freedom in a single-site block. 
(The number of degrees of freedom for the three- leg ladder is 2^ = 8, and that for 
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the effective Hamiltonian is 2^ = 4.) We also show the result of = 1/2 AFHC 
(effective Hamiltonian of the open ladder) or the open ladder for comparison. Note 
that S = 1/2 AFHC and the open ladders are expected to have power-low decaying 
correlation functions. 

First we show the result of the effective Hamiltonian (^. In Fig. |T] we check 
the truncation-dependent accuracy by calculating triplet excitation gap = 
-'-) ''total ~ 0) ^^^^ of truncation numbers m = 60 and 80 for the finite-size 

algorithm. By the figure we can see that m = 80 for the finite-size algorithm is 
enough for the calculation. It looks also that there is an excitation gap. The inset 
is the data for m = 80 and the fitting function 0.277 + 0.276 exp(-0.0387L). The 
extrapolated value of an excitation gap in the thermodynamic limit is 0.27(7). In 
Fig. ^ we plot the spin-spin and — correlation functions for the system size 
74 (m = 80). Both of them are exponentially decaying. The correlation length is 
estimated as 2.7(4). In Fig. |]we show the expectation value of the local Hamiltonian 
operator (|l3). The figure shows that the expectation value is not decaying. This is 
the strong evidence to show that the ground state of the system is dimerized and 
breaks the translational symmetry in the thermodynamic limit. To be convinced 
more carefully that all the system with < J2I J\ < 1/2, J^^oo'vs, in the same 
phase, we see the low energy spectra of the effective model with PBC by numerical 
exact diagonalization up to 12 sites (Fig. ^ shows 12 sites spectra for J2/ J\ = 0, 0.25 
and 0.5). There is no level crossing which supports no transition between J2/ Ji = 
and 0.5. Secondly we show the result of the three-leg ladder. We pick up the typical 
value of Jx = 1, 5 and 10 and the system size is 22 x 3. In Fig. |^ spin-spin correlation 
functions are plotted. They look like exponentially decaying. The correlation length 
tends to increase as J± decreases. We also check the truncation error by comparing 
the spin-spin correlation functions for truncation number m = 40, 60 and 80. (Inset 
of Fig. 1^) On the other hand, the data for the open ladders are power-law decaying 
as expected showing the system is critical. We show the data for the expectation 
value of the local Hamiltonian operator in Fig. These values do not decay, 
expressing long-range dimer-order. These results suggest that the cylindrical ladder 
is essentially represented by the strong coupling model (Hamiltonian(^) as expected 
by the renormalization group analysis. On the other hand, for the open ladders the 
expected results showing that the systems belong to the universality class of level-1 
SU{2) WZW model are obtained. Finally in Fig. |^we show the extrapolated values 
of the triplet excitation gap of the cylindlical ladder for typical values of J±/ Ji using 
data up to 22 X 3 sites. The fitting function is deduced from the effective Hamiltonian 
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(Fig. |T]), which is a + 6exp(— cL), where a, b, and c are fitting constants. 

§4. Higher Number of Odd Legs Cylindrical Ladders 

In this section we comment on cylindrical ladders with a higher odd number of 
legs in the strong coupling limit {J± oo). The ground states of an odd-sites 
antiferromagnetic ring are four-fold degenerate. (Though there is no exact proof 
of it, we expect this from the exact numerical diagonalization of small system size 
and field theory for a large number of system size. See ref. ^) We can derive 
the strong coupling effective Hamiltonian of the A^-leg [N = odd) cylindrical ladder 
(See Appendix for derivation of it.) which is given as, 

< = §E^ef^^'(«)- (11) 

i 

Here a is positive and dependent on the number of legs A^. For example 5-leg case 
a = 64/9. To calculate analytically the A^-dependence of a is too complicated for 
a higher number of legs. For small numbers of legs, however, we easily determine 
a numerically from the numerical exact diagonalization of system size N x 2. In 
Fig. ^ we plot the extrapolated value of a for A^ = 3, 5, 7, 9 and 11 using data up 
to J_L = 2048. The figure shows that a increases as the number of legs increases. In 
Fig. we plot the spin-spin correlation function and the local Hamiltonian operator 
with a corresponding to A^ = 5 and 7. The system size and truncation number are 
L = 50 and m = 80, respectively. These figures show that the effective models for 
cylindrical ladders of A^ = 5, 7 legs have dimerized ground states and it suggests 
that there is an excitation gap. 

§5. Summary and Discussion 

We study the ground state properties of the antiferromagnetic Heisenberg model 
on an odd number of legs cylindrical ladder. This model has frustrated interaction 
in the rung direction. The frustration induces the dimerization and the model has 
an excitation gap, that is different from the open ladder (which is with the open 
boundary condition in the rung direction.). The essential point distinguishing the 
cylindrical ladder from the open ladder seems to be the degeneracy of the ground 
state of odd-site antiferromagnetic Heisenberg chain i.e. the number of freedom of 
single site in the strong coupling effective Hamiltonian. First we consider the three- 
leg cylindrical ladder. We reduce the effective Hamiltonian (^ of the strong coupling 
limit (J_L oo), which is expected to be the fixed point of the renomalization group 
treatment. We introduce the next-nearest neighbor coupling (J2) for the effective 
Hamiltonian. For the model with special coupling (Ji = 2J2), we can construct the 
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exact ground state wave functions, which is direct products of two-site system wave 
function. This construction is closely analogous to Majumdar- Ghosh model. The 
wave function is apparently dimerized. The exactly solvable model is expected to 
have an excitation gap because the wave function is site-local. To confirm that the 
original system ( J2=0) has same nature with the exactly solvable model, we perform 
the numerical calculations using DMRG algorithm for the effective model and the 
three-leg cylindrical ladders. Adding to the spin-spin correlation function and triplet 
excitation gap, we calculate the expectation value of the local Hamiltonian operator 
which expresses the ground state being dimerized. The result shows that the three- 
leg cylindrical ladders are dimerized and essentially expressed by the strong coupling 
effective model. 

Lastly we briefly study the cylindrical ladders with a higher odd number of legs. 
We can calculate the strong coupling effective Hamiltonian for the general number of 
legs, which is the same form with the three-leg ladder. The A^-dependence is inserted 
into the parameter a. We have already known the exact solution for general a if we 
add the next-nearest neighbor interaction J2(= Ji/'^)- We expect the dimerization 
occurs for general even in the case of J2 = 0. We also confirm this expectation 
using DMRG calculations for small number of legs. We should comment on the 
ladder with a large odd number of legs. From the Fig. ^ we see that positive 
parameter a is monotonous as a function of A^ up to A^ = 11, we, therefore, expect 
that a for large number of legs are between 4 and 00. It is likely from Fig. ^that at 
least in the strong coupling limit the finite number of odd legs cylindrical ladder is 
dimerized in the ground state. The calculation in §4 is the strong coupling limiting 
case, therefore, in this study we cannot mention about the connection between 2- 
dimensional isotoropic Heisenberg model (which has Neel order and gapless spin 
wave excitation) and the effective model studied here. 
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Appendix: The Effective Hamiltonian for General 

In this Appendix, we derive the effective Hamiltonian (^) for general A^. As 
denoted in §4, the ground states of odd-sites antiferromagnetic ring are four-fold 
degenerate (two doublets) and these two doublets have opposite sign of momentum 
i.e. these momenta are k and —k. We write these four states as | j L), | | L), | | i?) 
and I i -R). Here L and R represent the momentum k and —k {k 0), and we set 
uj = exp{ik). Note that interchanging L and R corresponds to interchanging 00 and 
u^^. Matrix elements of S^^ and S'^j in the restricted Hilbert space are 

(T LI^^J T ^) = a„, R\Sl,\ ] R) = 

(T T R) = bn, (T RKJ ]L) = b: 



n' 



(T i L) = c„, (T I i R) = cl 

(TL|S+U/?)=rf„, {] R\St,\ I L) = d^, 

where n is the site number of odd site ring and a* represents the complex conjugate 
of a. It is obvious that interchanging | and | is interchanging the sign of a„ and 
bn and setting c„ = (i„ = 0. We set T is one-site shift operator for the n-direction, 
i.e. T""] ] L) = uj""] ] L) and T"| t ^) = cu"™] T R)- Then we conclude using the 
relation T'^S^^iT-'^ = ^ that 

a„ = ai, 6„ = 6ic.;-2("-i), c„ = Ci, rf„ = rfiu;-^^"-^). (A-1) 

The operator S'f of the effective Hamiltonian is the sum of 5*^^. Then using the 
relation eqs. (|A-l| ) we know a = 1/2 A. Setting the spin reversal operator P as 
n^=i2S'^j, we can derive the relations, 



ci = (T L\Sl,\ i L) = (i L\PSt,P\ ^L) = {i I T L) = (A2) 

i.e. Ci is also an real number. We also note that Y.n=i — ^n=i = Y.n=i '^n — 
Y,n=i dn = 0- Finally we obtain the effective Hamiltonian of the cylindrical ladder 
for general A, 

f ^ 

=Ji E ^S^S^i E K + l&nr(r+r- 1 + r,-r+ ,)) 

i I n=l 
^ n=l J 

In the derivation above, the SU{2) symmetry of the each interaction Sn,i ■ Sn,i+i 
is preserved, so it should be that the effective Hamiltonian (|A-3|) also possesses the 
symmetry. This means that Aaf = cf and |6ip/a^ = Idip/cf. Then we reached the 



result dn]) with a = |6ir/af (> 0). 
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We calculate a ior N — 5 as follows. The ground states for 5-sites ring {H — 
T,l=i Sn,i ■ Sn+i,i) with energy -\{2V5 + 3) are 

I t L) = \ (Alli) + /32|IIi)), 

where /?i = {u}~^ +a;~^)/2, /?2 = —{u}~^ +00) — | and ou — exp(|i7r). The vectors |Ii) 
and I III) are 

III) = ^ (I TTTU) + ^\ nut) + TUTT) + ^-'\ UTTT) + ^-'\ ilTTi)) 
liii) = (I Titit) + ^1 ititt) + oj'l Titti) + oj-'l ittit) + ^-'1 Ttiti)) ■ 

Other three states | i -^) , \ ^ R) and | J. R) are produced by interchanging f and J, 
or taking complex conjugate. After calculations we get, 

- 5(ifti^+ifen ■ (^■^> 

This leads to the result a = 64/9. 
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Figure captions 



Fig. 1. The finite size triplet excitation gap for the effective Hamiltonian {L = 74). The data 
is for the truncation number m = 60(0) and 80(x). The inset is the data for m = 80 and the 
fitting function 0.277 + 0.276 exp(-0.0387L). 



Fig. 2. The correlation functions (-l)''(S£/2'5i/2+r> (°) and (-!)'' (^1/2^1/2+^) (+) for the effec- 
tive Hamiltonian with system size L — 74 and truncation number m = 80 arc plotted in semi-log 
scale. The lines are for the guide to eyes (1.58x 10"^exp(-r/2.74) and 3.74x 10"^ exp(-r/2.74)). 



Fig. 3. The expectation value of the local Hamiltonian operator (— l)'^(ifioca'i^) (The symbol +, 
1 < r < L/2) is plotted in log- log scale for the effective Hamiltonian of the same L and m with 
Fig.2. The symbol O is for S =1/2 AFHC and the line is 0.20/ 



Fig. 4. The low-energy spectra for the effective Hamiltonian in the periodic boundary condition 
with J2/J1 = 0, 0.25 and 0.5. The system size L is 12. x-axis is lattice momentum in unit 2n/L. 
The symbols □, x, are spin-singlet with T^^^g^l = and r^otaj = 1 (The blacked symbol represents 
lowest two spin-singlet with T^^^^l = 0.), the symbols A, Q spin-triplet with r^ot^i ~ ^ and 
'^totai = 1' respectively. 
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Fig. 5. The correlation function (— l)'^('S'^/2'^L/2+r) three-leg ladders. (We use semi-log 

scale for cylindrical ladders and log-log scale for open ladders.) The system size is 22 x 3 and 
truncation number m = 60. We choose typical Jj_ as 1 (□), 5 (A) and 10 (x). The inset is 
the data for the cylindrical ladder with Jj^ = 5 with m = 40(x), 60(A), and 80(C>), we can see 
m = 60 is enough for the calculation. 



Fig. 6. The expectation value of local Hamiltonian operator (— l)'^(-ffioca'i ) (1 ^ ^ -^/2) for 
three-leg open ladders (black) and cylindrical ladders (white) with L = 22, N = 3. (The symbols 
O, □ and V a^re for Jj_ = 1, 5 and 10, respectively. The line 0.20/ ^/r are guide to eyes.) 



Fig. 7. The extrapolated value of the triplet excitation gap for the three-leg cylindrical ladders, 
using the data up to 22 x 3 with m = 60. The dotted line (0.27(7)) is the gap of the effective 
Hamiltonian (J_L — > oo). 



Fig. 8. The extrapolated value of a for the higher number of legs cylindrical ladder using the 
data of the numerical exact diagonalization of x 2 sites. 



Fig. 9. The correlation function {"^Y {S^/2'^L/2+r) '^^'^ local Hamiltonian operator 

{-ly {H[f^^i^) (1 < r < L/2) of the effective Hamiltonian with iV = 3 (0),5 (□) and 7 (A) 
are plotted in semi- log scale. 



15 



1.5 



T 1 1 1 1 1 1 1 r 

J^/J,=0 



o 



o 

□ 



1.0^1 



() 



0.50 



□ 



o 
□ 



0.0 



m- 




k (unit=27C/L) 



6 



A 



N 



N 



+ 



10" ? 



V 
I 

A 



10 ' ; 



+ 



10 ' : 



N 



N 

CD 
V 



2 10 ' . 



10 



-11 - 








-I I I I- 



10 



20 



30 



40 



50 



r 



10° 



A 



N 

(/5 



N 

V 



CM 



10 ^ - 



10 




10 



10" 



cylindrical ladders 



10' 



10^ 



□ □ 



□ □ 



10' 



10* 



10 



12 



r 



0.30 




0.00 







6 

J 



8 



10 12 



20.0 



15.0 - 



10.0 - 



5.0 - 



0.0 








6 8 10 12 



T 1 1 r 



T 1 1 r 



T 1 1 r 



T 1 1 r 



A 



^aaaaaaaaaaaaaaaaaaaaaaaa 

^□□□□□□□□□□□□□□□□□□□□□□□n 



.^OOOOOOOOOOOOOOOOOOOOOOOO 



10 ' I 




(-1) <H''^' > 



local 



10 ' . 



10' . 



o 



(-1)'<S^ > 

L/2 L/2+r 



10" . 



o A 

R 







10 ' . 







-I I I u 



10 



15 



-I I I I u 



20 



25 



30 



r 



1.5 



T 1 1 1 1 1 1 1 1 1 1 r 

JJJ^= 0.25 



1.0 - 



0.50 



o 
□ 



o 
□ 



o 
□ 



o 



□ 



[] 



0.0 



m- 




k (unit=27C/L) 



6 



1.5 



T 1 1 1 1 1 1 1 1 1 1 r 

J^/J,=0.5 



1.0 - 



[3 



□ 



O 



□ 



O 



□ 



o 



□ 



0.50 



0.0 



m- 




k (unit=27C/L) 



6 



